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Abstract 

We study nonlinear dynamics in models of Lorentz-violating massive grav- 
ity. The Boulware-Deser instability restricts severely the class of acceptable 
theories. We identify a model that is stable. It exhibits the following bizarre 
but interesting property: there are only two massive propagating degrees of 
freedom in the spectrum, and yet long-range instantaneous interactions are 
present in the theory. We discuss this property on a simpler example of a pho- 
ton with a Lorentz-violating mass term where the issues of (a)causality are 
easier to understand. Depending on the values of the mass parameter these 
models can either be excluded, or become phenomenologically interesting. We 
discuss a similar example with more degrees of freedom, as well as a model 
without the long-range instantaneous interactions. 



1 Introduction and summary 



Models of a massive/metastable graviton could shed some light on the cosmo logical 
constant problem (see, e.g., jTj). Given the ultraviolet problems of gravity, provi- 
sionally, one would like to find a model that could be regarded as a classically and 
quantum mechanically consistent low-energy effective field theory (for a description 
of gravity as an effective theory see Ref. j2]). In searching for such a theory, that 
could also preserve Lorentz invariance in four-dimensions, one typically encounters 
the following three major problems: 

Problem 1. Linear discontinuity. Known as the van Dam-Veltman-Zakharov 
(vDVZ) discontinuity • To understand this problem let us ignore for the time 
being all possible nonlinear self-interactions of a spin-2 field. The Lagrangian of this 
theory without ghosts and tachyons was uniquely determined by Fierz and Pauli 
(FP) [Sj. Irrespective of the details of the Lagrangian, Lorentz invariance dictates 
that a massive spin-2 state has to have five physical polarizations, as opposed to 
a massless graviton with only two polarizations. This is what gives rise to the 
vDVZ discontinuity: one out of the extra three degrees of freedom couples to the 
trace of the stress-tensor, and, no matter how small the graviton mass, gives rise to 
experimentally unacceptable predictions either for the light bending by the Sun or 
for Newtonian interactions. Thus, the FP Lagrangian gives a theoretically consistent 
model of a massive spin-2 state without self-interactions, however, if this spin-2 state 
is declared to be a graviton, the model is in contradiction with observations. 

Problem 2. Strong coupling. The above arguments were based on the linearized 
approximation. In a theory of gravity we should take into account nonlinearities. It 
was first observed by Vainshtein ,6j that the arguments leading to the vDVZ discon- 
tinuity fail once the nonlinear interactions are taken into account. This is because 
linearized approximation breaks down, and, to make predictions within the Solar 
system one should either solve the nonlinear equations exactly, or come up with an 
alternative viable approximation. The latter two approaches can restore the con- 
sistency of the predictions of classical massive gravity with observations [HIIZI. The 
breakdown of the linearized approximation takes place because of the nonlinear but 
classical self-interactions of the extra polarizations |7j. This suggests a problem for 
quantum theory where the same nonlinear interactions appear in the loop diagrams, 
leading eventually to a very low ultraviolet cutoff |H| for a quantum graviton scat- 
tering on an empty background. Moreover, quantum loops are expected to generate 
higher dimensional operators that are suppressed only by this observationally un- 
acceptable low cutoff Note that in a theory of gravity a cutoff has no universal 
meaning. For instance, for graviton scattering in a background of classical sources 
the physical cutoff would depend on local curvature invariants. This in principle 
could be used to try to overcome the strong coupling problem (along the lines of 
jni E]) if the model were consistent otherwise. This brings us to the issue of a 
nonlinear instability of the FP model. 
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Problem 3. Nonlinear instability. Also known as the Boulware-Deser (BD) in- 
stability. From our standpoint, clarified below, this is the most severe problem. 
It emerges already at the classical level. To quickly sketch the essence of the 
BD instability one can look at a scalar field model with the Hamiltonian H = 
{dt(j)Y + iy(j)Y + 'm'^4>^ + cr(((9t0)^ + (V0)^) + m^cr^0/2. In the quadratic approxi- 
mation this describes a single free scalar field (f). However, because of the nonlinear 
(cubic in this case) terms the instabilities set in. One can simply integrate out the a 
field and obtain that the Hamiltonian contains the term — + (V</')^)^/2m^(/), 
which is sign-indefinite and unbounded below. This shows the essence of the BD 
instability in an oversimplified model (for complete treatment, see section 3). In this 
toy model the problem can be cured by adding new terms into the Hamiltonian. One 
might hope for a similar remedy in the FP theory. Indeed, a nonlinear completion of 
the FP gravity is not unique and, in any event, one should expect quantum loops to 
generate all sorts of new nonlinear terms. There are examples of nonlinear systems 
where certain classical instabilities are removed by quantum-loop-generated terms^. 
However, to the best of our knowledge, there is no evidence for this to be happening 
in the FP model. In particular: (a) The worst part of the BD instability exists 
because of the nonlinear interactions between the tensor sector of the conventional 
Einstein-Hilbert Lagrangian and the extra polarizations of a massive graviton that 
are only present in the mass term (i.e., the Nambu-Goldstone (NG) sector of massive 
gravity), (b) It was shown in Refs. P^ IT^ that the BD type instabilities persist for 
an arbitrary polynomial in fields completion of the FP gravity. Note that these terms 
include all sorts of derivatives of the Nambu-Goldstone boson, (c) One might hope 
that these instabilities go away once the terms that contain derivatives of fields are 
included. However, this would require an infinite number of fine tunings for which 
no symmetry principle is known in the four- dimensional context ^. The unstable 
solutions found in Ref. [T] were of a cosmological type. It would be interesting to 
conduct similar studies for spatially inhomogeneous localized sources. 

Summarizing, the BD instability is the most severe problem: unlike Problems 
1 and 2, it questions the very consistency of a theoretical model itself. So far no 
concrete cure was proposed in the context of a Lorentz-invariant local field theory 
with a finite number of degrees of freedom and without an infinite number of fine 
tunings. Surely, a consistent theory of nonlinearly interacting massive spin-2 should 
exist, it just has not been formulated yet. Therefore, until shown otherwise, we will 
assume that the models with the DB instability should be avoided. 

In an ideal case, one would like to have a model in which all the above three 
problems are absent. However, if one should compromise between Problems 1,2 and 

-'^For instance, some models exhibit chaotic behavior at the classical level while the chaos is 
eliminated by quantum corrections. 

^One can also try to modify the linear part of the FP theory by introducing heavy states at the 
UV cutoff of the theory. Explicit calculations, similar to those of section 3, show that the rapid 
BD instability persists in this case too. We thank Nima Arkani-Hamed for suggesting to study 
this question. 
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3, as it will be the case in one particular examples below, our approach will be to 
worry first of all about the Problem 3. This is because Problems 1 and 2, although 
unpleasant from the point of view of practical calculations, can be taken care of 
consistently. For instance, in the DGP model ^2] these problems are solved at the 
classical [7^ as well as at the quantum level P ITU]^. 

Recently, a new approach to massive gravity was initiated in Ref. J7j: the idea 
is to give up Lorentz-invariance which could be spontaneously broken when graviton 
acquires its mass |T7j. Subsequently, in Ref. 18 a general parametrization of the 
Lorentz- violating (LV) graviton mass term was proposed and the models evading the 
Problems 1 and 2 were identified. More general studies of the proposal of ^H] were 
performed in Ref. The discussions in Ref. [18 were restricted to the linearized 
theory. The purpose of the present work is to study a complete nonlinear dynamics 
in the models with the LV mass terms, and in particular address the Problem 3. We 
will find that many of the LV mass models suffer from the BD instability. However, 
there are at least two classes of models that can evade the problem. The first 
class exhibit surprising properties: even though there are only massive propagating 
degrees of freedom, the models exhibits a long-range instantaneous interactions. 
This could be phenomenologically deadly or interesting depending on the value of 
the graviton mass. It is very likely that the properties of some of these models won't 
be affected by radiative corrections since they are protected by certain symmetries. 
(A different model but with somewhat similar properties was discussed in jTH] in the 
linearized approximation, see also (2iij). The second class of the BD stable models 
contains all the massive degrees of freedom and no long range interactions. However, 
the issue of radiative corrections for these models remain open, without an infinite 
number of fine tunings, these models are likely to exhibit the instabilities at the 
quantum level. 

The work is organized as follows. In section 2, as an instructive example, we 
discuss a Lorentz-violating theory of a massive photon. This model contains only 
massive propagating degrees of freedom (two massive polarizations of an electro- 
magnetic wave), and nevertheless, there are long-range instantaneous interactions 
in the theory. We briefly discuss whether this type of models can be consistent with 
observations. In section 3 we overview the BD instability in the FP theory and show 
how it also appears in the LV models. In section 4, first we discuss a model that has 
no BD instability and propagates only two massive degrees of freedom (transverse 
polarizations of gravitational waves). We study the long-range interactions in this 
model, showing that there is no vDVZ discontinuity. Finally we discuss two other 
models, one with 5 massive degrees of freedom and the long-range-interactions, and 
another one with six degrees of freedom where all the interactions are screened. 

^We also note that there exist models that have no Problems 2 or 3 |14[ I15| . however they 
exhibit the vDVZ discontinuity. A version of 13 discussed in 16 could potentially evade all the 
problems in the weak coupling regime, however some nonlinear issues should still be understood 
in that approach (see the discussions section in |lf)j\ 
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2 Warming up with photons 



As a toy but very interesting example we consider QED with Lorentz- violating mass 
term for a photon^ (for convenience we call it QED') 

/ = - \m^A,A^ - A,J^ , (1) 

where /i, = 0, 1, 2, 3; i,j = l, 2, 3; is a conserved current 9^ = 0, and the mass 
term breaks explicitly the Lorentz group down to the group of three dimensional 
spatial rotations 5*0(3) (our choice of the Lorentzian signature is [ — h ++]•)• One 
could think of this model as arising from a Lorentz-invariant theory in which certain 
fields acquire Lorentz- violating VEV's (see, e.g., [T7j, [22] )■ These VEV's set a 
preferred frame in the Universe in which the model (QJ is defined. 

The Lagrangian is invariant under spatially independent gauge transforma- 
tions of the fields 6Ao = do\(t), 6Aj = 0. Because of this and conservation of 
no new terms are generated by quantum loops in (^. The equations of motion of 
the model are 

d''F,^-m''5;A = J^. (2) 

There are two key properties that follow from First, the zeros component of 
this equation implies that the Gauss's law is identical to that of QED: 

d^E, = -Jo, (3) 

where Ej = Fqj is an electric field and Jq is a charge density. Hence, in QED', like 
in QED, the electric field is not screened). Second, taking the partial derivative of 
both sides of Q we find a Lorentz-violating analog of the Proca condition 

d'A = 0. (4) 

As a result, there remain only two dynamical degrees of freedom in the theory: Aq 
is not dynamical and one of the A/s can be expressed though the other two using 
(jD). Both propagating degrees of freedom are massive. 

However, this is not all. We note that (jH) coincides with the Coulomb gauge 
fixing condition of QED. Because of this, a free photon propagator of 

^oo(A:) = ^, A,(fc) = f% - ^] , . , (5) 

resembles a causal Coulomb gauge QED propagator. The physics of QED' ((H) is 
rather different, however. Like in QED, there is an instantaneous Coulomb potential 
(the zero- zero component of the propagator). This component has no imaginary 

^Michele Papucci and Matthew Schwartz also studied this model for a photon. We thank 
Matthew Schwartz for communications on this. 



5 



part, hence, there are no physical degrees of freedom mediating the instantaneous 
potential. The spatial components, on the other hand, have an imaginary part. 
This corresponds to the two physical degrees of freedom. Unlike in QED, in the 
present case both of these degrees of freedom are massive. This has a dramatic 
consequence: as we will see shortly, the instantaneous potential is not canceled in 
physical observables for time dependent sources. The remaining instantaneous field 
is small for transverse sources with a typical momentum/frequency scales 3> m, but 
can become essential for scales ^ m generating action- at- a- distance. To examine 
this question in detail we follow closely massless QED in Coulomb gauge. In the 
latter case one postulates (jH) as a gauge condition. As a result of this 

Mr,t)^± fd\''j^, (6) 

4tt J \r — r\ 

is an instantaneous potential. The expression ® is identical in QED and QED'. On 
the other hand, the equation for the vector potential differ in the two models. The 
spatial part of (0) reads: 

{d^ - m^)A^ = J, - d.doAo = Jf . (7) 

The mass term on the l.h.s. is present only in QED' but not in QED. Both in QED 
and QED' the vector potential Aj has an instantaneous parts. In QED this part 
exactly cancels © in physical observables such as the electric field Ej = —OjAq + 
doAj. However, this cancellation is not exact in QED'. To see this we write: 

A,(r, t)= I dh'dt'DRir -r'-t~ t')Jf{r\ t') , (8) 
where the retarded Green's function 



is expressed through the step function 6, Dirac delta function 6, and Bessel function 
J' Note that the mass dependence in the Green's functions is additive. We 

used this to denote the massless function by D^^^ and the addition due to the mass 
term by D^. Using ^ we can write 

A,(r, t) = Af^^'ir, t) + J d^r'dt'D'^ir -r';t- t')jf{r\ t') , (10) 

where A^^^{r,t) is a vector potential of massless QED in Coulomb gauge. The 
latter, as we mentioned before, cancels the instantaneous scalar potential © to 
produce the retarded electric field E^^^. Therefore, 

Ej{r, t) = Ef^^ir, t) + do J d^r'dt'D'^^ir -r';t- t') Jf (r', t') , (11) 
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and the instantaneous part is now contained only in the last term of this expression. 
The latter can be calculated as follows: 



I 



I 



rf^fce^^" Re (e*^'^'+'"'* Jf{k, ^fW+m^) - e'^''^' Jf{k, \k\)) 



(12) 



where 




(13) 



is the Fourier transform of the transverse current^. 

In general, the expression (fT^ is nonzero even for t = 0. It appears that an 
information from an event that took place at t = 0,r = can instantaneously be 
transmitted to a point that is far away from this location. This gives rise to the 
action-at-a-distance. The question how important this instantaneous interaction 
is depends on properties of sources. For a transverse source of a typical momen- 
tum fco and typical frequency ujq the effect is negligible as long as kQ ^ m. For 
ko m the effects can be appreciable when tuo ~ m or cuq ~ fco- In this case the 
instantaneous electric field would decay with distance as ~ 1/r. If the source has 
no typical frequency, juts a typical momentum, then the instantaneous interactions 
will be important for ko <^ m, and vice versa, for a source of a typical frequency c^o 
and no typical momentum the dangerous interactions will be present for ujq ^ m. 
In practice, to produce a low-momentum/frequency signal that could trigger the 
instantaneous interaction, will itself take certain characteristic time. It would be 
interesting to study the phenomenology of these interactions for realistic sources to 
put bounds on m [23] ■ Note also that magnetic field has no instantaneous parts in 
QED': B = curM and the curl eliminates the instantaneous part of A. 

The presence of the long-range interactions can also be understood from the 
Hamiltonian formulation of (Q) where Aq acts as a Lagrange multiplier: 



Here Pj = Fjq is the canonical momentum. Variation w.r.t. Aq gives rise to the 
Gauss's law, d^Pj = Jo, which is identical to the Gauss's law of massless QED. This 
guaranties that the theory possesses the long-range interactions in spite of the fact 
that the two dynamical propagating degrees of freedom are massive^. 

To summarize, there are two massive propagating degrees of freedom, however, 
there still exists a long-range instantaneous interaction in the theory. Depending 
on the value of m this can exclude a given model, or be potentially interesting for 
phenomenological applications. 

^For simplicity (|12(l is written for a source such that J{k,—u!) = J(k,uj). However, a general 
expression can readily be obtained. 

^The condition Aj = is obtained by taking a derivative of one of the Hamiltonian equations 
of motion. 



^ = 1^! + k^vkdjAkf + A] + A, {d,P, - Jo) + A, J, . (14) 
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One can of course modify the model (P) by adding a "mass term" a m'^Al with 
some nonzero positive coefficient a. In this case the long-range interactions are 
removed from the theory (a = 1 corresponding to a massive Lorentz invariant pho- 
ton). However, it is worth pointing out that a = is an enhanced gauge symmetry 
point and unlike the a 7^ cases should be stable under radiative corrections^. 

3 Nonlinear instabilities in massive gravity 

In this section we summarize how nonlinear instabilities appear in a Lorentz-invariant 
theory of massive gravity (the FP gravity) JT] , and show that the similar instabili- 
ties exist in many of the LV mass models once the nonlinear interactions are taken 
into account. We will also find the conditions under which these instabilities can 
be removed in a Lorentz-violating massive theory. In the latter case, however, one 
typically ends up with long-range interactions, similar to those studied in the previ- 
ous section. In this respect the theory is half-massive: all its degrees of freedom are 
massive nevertheless there are long-range interactions. One atypical example that 
evades the instabilities and long-range instantaneous interactions will be given at 
the end of the next section. 

We first review briefiy the Hamiltonian construction of JT] to identify the terms 
that are responsible for the instabilities. Then, we will remove these terms in a 
Lorentz-violating theory. Let us start by a brief reminder of the ADM formalism 
This would be a natural formalism for nonlinear formulation of the LV mass 
gravity. Consider a foliation of space-time by hyper-surfaces parametrized by 
a time variable t. The four-dimensional metric is replaced by the following three- 
dimensional variables: 

7., = ^?.,; iv = (- (V°)-^/^; N,^ (^W. (15) 

In terms of these variables the invariant interval takes the form: 

ds'^ = -{N^ - NjN^)dt^ + 2Njdx^ dt + -Hjdx' dx^ , (16) 

where all the spatial indices are contracted by means of the three-dimensional metric 
7jj. is called the lapse function and Ni is the shift function. With these definitions 

^- Wg = N^, (^)/?= ^'^R + K,,K'^ -K\ (17) 
where Kij is the extrinsic curvature of St 

K,, = ^AT-i [7,,. - V.AT, - V,N,] , (18) 

^Thc a — 1 case is also stable because of the restored Lorentz invariance and conservation of 
the current. 
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and K is its trace. The extrinsic curvature is related to the canonical momentum 

61 



TC 



^{K'^-KY^). (19) 



■ 6iij 

The Hamiltonian of the Einstein gravity in terms of the above variables reads 

Heh = n^'iij -l = Vl [nR' + NiR"] , (20) 

where 

i?o = - + 7-i(7ri,7r'^' - ivr^) , K = -2V,{^-'^^7r'') . (21) 

and Ni appear linearly in the Hamiltonian ((201). Hence they are Lagrange multi- 
pliers variation w.r.t. which gives the constraints R^ = and R^ = 0; Heh = on 
the surface of the constraints. 

Let us now turn to massive FP gravity for which the mass terms is written as 

- Wihl - hf) = -Im' \hl -h'- 2Nf + 2hil -N'- N')] , (22) 



where the second equality is obtained by expressing h^^, = g^^ — r]^y in terms of 7jj, 
N and Ni (note that hij = '-/ij — rjij, and h = Y^'hij). The key role is played by the 
terms in ()22|) which are quadratic in and Nj. Because of these terms N and Nj 
ceases to be Lagrange multipliers in the massive theory. Variations w.r.t. and Nj 
lead to the following equations: 

^7 = 2m^hN , = 2m\r]'^ - hY^)Nj . (23) 

These are not constraint equations any more but serve to determine and Nj. 
Substituting these solutions into the full Hamiltonian we obtain: 



H = ^ 1^^^^^ + iR'W' - hYr'R'] + -h' + 2h) 



(24) 



This is a Hamiltonian of an ill defined theory. The first term on the r.h.s. is 
unbounded below and singular in m and h. For instance, consider ^/jR^ fixed and 
i?* = 0; when h ^ 0^ the term in the Hamiltonian density Ti. ~ {^R^Y / {m?h) 
is not bounded below. This demonstrates the presence of a ghost-like instability in 
the theory. This instability can manifest in many ways even at the classical level, 
and the time scale of the instability can be very short jl2j. Such a theory is hard to 
make sense of. 

The BD problem is associated with the terms that in the linearized theory looks 
as /loo^j- The mass term of the model analyzed in Ref. ^Hl is 

Lm = ^ (mlhlo + 2mlhlj - m^/i^. + ^2^2 _ 2mlhooh) , (25) 
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with niQ = and all the other mass parameters being nonzero with a certain hierar- 
chy between them ^Hj. A straightforward non-linear completion of the mass term 
gives rise to the BD instability. This is because /iqo = 1 ~ + and N 
ceases to be a Lagrange multiplier, as in the FP gravity. While in the FP model this 
instability cannot be removed in a Lorentz-invariant fashion, in the present context 
Lorentz-invariance is broken anyway, and nothing prevents us to eliminate the dan- 
gerous term by judicially choosing to be zero. This choice is a point of enhanced 
gauge symmetry and should be stable under loop corrections. However, the physics 
of the model with mo = m4 = is dramatically different - there appear long-range 
interactions. This will be discussed in detail in the next section. 



4 Stable models 

1. Half-massive gravity 

We consider a simple Lorentz-violating generalization of the FP mass term: 

Ah = -\m^^N{hl-ah''), (26) 

where a is a constant. As before, all the indices are contracted by 7,^-. We think of 
this theory as being obtained from a Lorentz-invariant model through spontaneous 
generation of a preferred frame, similar in spirit to |17j . 
The total Hamiltonian in this case takes the form: 



The constraint equations that follow are: 



N[B!^ + ^m\hl - all') ) + NiR 



(27) 



i?o = --m^{hl - ah^) , R' = 0. (28) 
2 

The Hamiltonian ()27p on the surface of constraints ()28|) is zero, just like in the 
Einstein theory, hence, the BD instability is gone. The model is invariant in 
the linearized approximation under coordinate-independent gauge transformations: 
6h^i, = d^C,y{t) +d^(^(t), as well as under the transformations with a gauge function 
= {Co{t, = 0) (at the nonlinear level there is a symmetry w.r.t. the spatially 
independent transformations). Due to this, we expect that the properties of this 
model will stay stable under quantum loop corrections. 

Let us now couple this theory to matter. The {00} component of the equation 
of motion takes the form (we use the units Mpi = 1): 

+ \m\hl-ah') = Too. (29) 

The {Oj} equations are identical to those of the Einstein theory and read as Rj = 
2Toj. Finally, the {ij} equations are: 

Gij + ^m^ihij - a6ijh) = Tij . (30) 
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We see that in the hnearized theory the {00} equation is also identical to 
the linearized massless Einstein equation. Therefore, it is only the {ij} equation 
that differentiates (|26p from the Einstein gravity in the linearized approximation. 
Because of this one should expect the vDVZ discontinuity to be absent. This can 
be checked in a rigorous way. Let us follow the decomposition of Ref. 



/loo = i', (31) 
hoi = Ui + diV , (32) 
hij = Xij + d(iSj) + didja + 6ijT , (33) 

where Xtj is a transverse-traceless tensor, sj is a transverse vector while the other 
fields are scalars. The gauge invariant combinations are: a tensor Xij^ a vector 
Wi = Ui — doSi, and two scalars r and ^ = ip — 2dov + d^a. The conventional 
coupling to a conserved matter stress-tensor h^^T^'^ can be written in terms of these 
invariants: 

XijTij - ^WjToj + $Too + rTjj . (34) 
Solving the corresponding linearized equations for a 7^ 1 we find: 

—Oq + a — 

^ = 2^'^°°' ^ 1"^°^'' ^^^^ 

(Here T^j stands for a transverse-traceless part of the tensor). These expressions 
give exactly the fields of the Einstein theory in the limit m —>■ 0. Therefore, there is 
no vDVZ discontinuity. Note that the only propagating degrees of freedom are two 
polarizations of the transverse-traceless tensor Xij- The spectrum is free of ghost 
and tachyons. Similar properties have been found previously in a stimulating work 
|20j . where a somewhat different model was discussed^. 

On the other hand, the above system exhibits the same type of instantaneous 
interactions as QED' discussed in Section 2. This is because the instantaneous parts 
cannot be exactly canceled as long as there is a mass term in the denominator of 
(jSni)- However, as we discussed in the gauge field case, these instantaneous inter- 
actions can only be probed by very low- momentum/frequency sources. It would be 
interesting to study the phenomenology of this model 



2. More degrees of freedom 
A generalization of the above model can be obtained by adding to (j26p a mass term 



^However, non- linear stability of the model of Ref. (201, which is a suspect in the context of the 
discussions of our section 3, has not been studied in Ref. |20| . 
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for Nf 



Ak = cm'^Nt, (38) 

where c is a constant. Doing so we add three additional degrees of freedom to the 
theory. The Hamiltonian now takes the form: 



N (r^ + -m^{hl - ah^)] + Ni (i?^ - m^cN 



(39) 



2 

The corresponding constraint equations are: 

1 

'2' 

R = 2m'^cN\ (41) 



i?o = --m\hl-ah^), (40) 



As long as c 7^ the only true constraint is (jlOl), since PT|) does not restrict the 
number of propagating degrees of freedom but acts as an algebraic equation deter- 
mining A^-^ . Therefore, this counting tells us that the number of propagating degrees 
of freedom is five. Solving ()41|) for A^* and substituting the result in 7i we find the 
Hamiltonian 

^ = v^Z^ ' (42) 
4m^ c 

which is positive semi definite as long as c > 0. This model also has no vDVZ dis- 
continuity. The calculations are similar to those presented above but more tedious, 
the spectrum contains no ghost of tachyons. The massless limit of ()39|) is regular 
and one recovers in this limit the Einstein gravity. Moreover, the theory is symmet- 
ric w.r.t. spatially independent transformations of the time variable, and, exhibits 
the instantaneous interactions. Further interesting properties of this model will be 
discussed in 



One can also evade the BD instability and the presence of the long range inter- 
actions by adding into the Lagrangian yet another term 

Ah = m^^P^iN), (43) 

where P2{N) is a polynomial in A^ of degree 2, namely P2{N) = cq + CiN + C2A^^. 
We require that there are no constant and linear terms in the linearized Lagrangian. 
This gives the relations cq + ci + C2 = and C2 -|- |ci = 0, with a solution ci = — 2co 
and C2 = Cq. Hence, P2{N) = Cq{N — 1)^. The Hamiltonians of the system is that 
of the previous examples plus the new term AH = —rn?^P2{N). The resulting 
constraint equations are: 

R"" + Kn\hl-ah^) = m^P^{N), (44) 
R' = 2m^cN\ (45) 
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Solving the constraints w.r.t. and A^* we obtain 



N 



2im?c 



and the Hamiltonian 



where Rl^^^ 



dO _|_ 
-'"'mod 



(-^mod)^ _|_ 



(46) 
(47) 

(48) 



R^ + ^rn?{h1j — ah?). The above Hamiltonian is not necessarily 
positive semi definite, however, it is bounded from below as long as c, Cq > 0. There 
are six degrees of freedom propagating in this model. This could exclude the model 
based on the Solar system data. However, it is not impossible to imagine that 
the nonlinear effects suppress the couplings of the extra polarizations to matter at 
observable distances, in analogy with [71I2S]. A more serious problem of the model 
(jIHj) is the absence of a symmetry principle that would guarantee the stability w.r.t. 
quantum corrections. 
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